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THE  EFFECTS  ON  STRESS  SINGULARITIES  OF  INTRODUCING  COHESIVE 
STRESS-SEPARATION  LAWS  AS  BOUNDARY  CONDITIONS  FOR  ELASTIC 
PLATES  IN  EXTENSION 


Introduction 

The  stress  fields  in  angular  elastic  plates  in  extension  were  first  treated  by  Knein 
[1],  and  subsequently  systematically  identified  by  WiUiams  [2].  These  studies  reveal  the 
possibility  of  elastic  stress  singularities  at  the  vertex  of  the  wedge.  For  example,  a  stress- 
free  wedge  of  angle  2k  realizes  a  model  of  a  crack  having  stresses  which  behave  as  the 
inverse  of  the  square  root  of  the  distance  from  the  crack  tip.  Such  fields  are  at  variance 
with  file  assumptions  underlying  the  theory  from  which  tiiey  came.  Consequently,  care 
needs  to  be  exercised  in  interpreting  them.  Even  with  care,  resulting  inferences  with 
respect  to  structural  integrity  are  not  always  as  reliable  as  desired  in  practice  (see,  for 
example,  [3]).  Accordingly  it  is  of  some  interest  to  explore  alternative  models  which  are 
free  of  these  nonphysical  singular  stresses.  This  report  considers  a  means  to  this  end, 
namely  the  introduction  of  cohesive  stress-separation  laws. 

Barenblatt  was  first  to  introduce  cohesive  stress-separation  laws  to  rid  cracks  of  their 
stress  singularities;  an  extensive  account  of  his  work  is  given  in  [4].  Essentially  his 
approach  cancels  the  opening  singularity  produced  by  loading  remote  from  the  crack 
with  the  closing  singularity  produced  by  cohesive  stresses  on  the  crack  flanks  near  the 
crack  tip. 

To  explain  further,  we  consider  the  example  of  a  central  crack  in  an  infinite 
elastic  plate  under  far-field  uniform  tension  (Fig.  1).  The  crack  has  length  2a,  the 
applied  tension  magmtude  Oq.  To  describe  the  configuration,  we  use  rectangular 
Cartesian  coordinates  x,  y,  with  origin  O  at  the  rightmost  crack  tip.  Then  the  normal 
definition  of  the  stress  intensity  factor,  Ki,  associated  with  this  crack  tip,  has 


Kj  =  lim  Oy  I  , 

y=0 


(1) 


2 


Fig.  1.  Crack  with  local  cohesive  zone  (after  Barenblatt  [4]) 
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where  Oy  is  the  normal  stress  component  in  the  y-direction.  Accordingly,  Kj  is  the 
coefficient  of  the  aforementioned  inverse-square-root  singularity  in  the  crack-tip  stress 
field.  Its  values  for  the  central  crack  configuration  of  Fig.  1  are  as  follows.  For  the 
applied  stress  Cq,  we  have  (see,  e.g.,  Tada,  Paris  and  Irwin  [5],  p.  5.1) 


Ki  =  OoVita . 


(2) 


For  the  cohesive  stress,  Oc  =  ac(4)  where  ^  =  -x,  we  have  (see,  e.g.,  Tada  et  al.,  [5],  p. 
5.11) 


U  0  T/«2a  -4)' 


(3) 


where  Aa  is  the  extent  of  the  zones  in  which  the  cohesive  stresses  act.  Thus  by  ensuring 
the  magnitudes  of  the  stress  intensity  factors  in  (2),  (3)  are  equal,  the  advertised 
cancellation  takes  place  and  effects  a  singularity-free  crack. 

In  addition  to  introducing  cohesive  stresses,  Barenblatt  made  some  assumptions 
regarding  their  distribution.  First,  he  assumes  that  the  extent  of  the  cohesive  zones  is 

small  ([4],  p.  78),  i.e.,  Aa  /  a  « 1.  Under  this  assumption,  the  stress  intensity  factor  of 
(3)  is  given  by 


(4) 
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as  Aa/a—^0,  wherein  O  is  the  large  order  O  symbol.  Barenblatt  terms  K  of  (4)  the 
modulus  of  cohesion.  His  second  assumption  then  has,  in  effect,  that  the  maximum 
possible  value  of  K  does  not  depend  upon  the  applied  loading,  Oq,  and  is  always  the 
same  for  a  given  material  ([4],  p.  79).  Thus  when  this  maximum  is  attained  and  fracture 
is  about  to  occur,  provided  the  configuration  is  maintained  singularity  free,  the  stress 
intensity  factor  due  to  Oq  has  also  attained  a  critical  value.  Essentially,  then,  tiie  stress 
intensity  factor  due  to  applied  loading  becomes  the  key  parameter  controlling  fracture, 
the  same  fracture  criterion  as  used  in  fracture  mechanics  when  singularities  are  present. 
Goodier  [6]  and  Rice  [7]  also  provide  arguments  that,  as  a  result  of  his  assumptions, 
Barenblatfs  approach  reduces  to  the  same  as  for  cracks  with  singularities  present. 
Hence,  while  Barenblatt  does  indeed  cancel  singularities  for  cracks  by  introducing  the 
concept  of  cohesive  crack-flank  stresses,  the  manner  in  which  he  does  so  leads  to  an 
approach  which  is  equivalent  to  that  used  when  singularities  are  present  as  far  as 
fracture  goes.  As  a  consequence,  this  cancellation  of  singular  behavior  is  perhaps  not 
that  significant  from  a  practical  viewpoint. 

Here,  in  contrast,  we  do  not  assume  the  cohesive  zone  to  be  necessarily  small.  In 
fact,  we  let  it  extend  ahead  of  the  crack  tip,  as  indicated  in  the  sketch  of  Fig.  2.  This 
extension  changes  the  boundary  conditions  ahead  of  the  crack.  Since  boundary 
conditions  play  a  critical  role  in  determining  the  possible  singular  behavior  present,  this 
change  raises  the  question,  "what  type  of  singularity  can  the  crack  now  have?"  We  seek 
to  investigate  this  question  herein,  both  for  the  cracked  plate  in  particular,  and  for  the 
angular  plate  in  general. 

We  begin  in  Section  1  with  a  formal  problem  statement.  Next,  in  Section  2,  we 
develop  boundary  conditions  equivalent  to  cohesive  stress-separation  laws  which,  in 
turn,  lead  to  the  eigenvalue  equations  given  in  Section  3.  We  then  examine  some 
associated  eigenfunctions  for  special  cases  in  Section  4,  and  close  by  offering  some 
concluding  remarks. 
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1.  Formulation 

Here  we  specify,  in  some  detail,  the  class  of  angular  elastic  plate  configurations 
to  be  examined  for  possible  stress  singularities. 

To  facilitate  the  treatment  of  the  angular  plate  geometry  (Fig.  3),  we  use  polar 
coordinates  r,  0,  sharing  the  same  origin,  O,  as  their  rectangular  coimterparts  x,y,  and 
related  to  them  as  in: 

X  =  r  cos  0,  y  =  r  sin  0,  (5) 

for  0<r<oo,  O<0<  Itc.  Then  the  region  of  interest,  SR,  is  given  by: 

SR  =  {(r,0;|O<r<oo,  O<0<<1)},  (6) 

where  <1)  is  the  vertex  angle  of  the  angular  plate.  With  these  geometric  preliminaries  in 
place,  we  can  formulate  the  class  of  problems  of  concern  as  follows. 

In  general,  we  seek  the  two-dimensional,  elastic,  stress  components  0^,09, x^e, 
and  associated  displacements  Uj.,Ue,  satisfying  the  following:  the  stress  equations  of 
equilibrium  in  the  absence  of  body  forces, 

aor  1  dXre  Or-ae 
dr  r  dQ  r 

(7) 

1  dcQ  3Xre  ,  2Xre  _ 

r  30  ar  r 

on  SR;  the  stress-displacement  relations  for  a  homogeneous  and  isotropic,  linear  elastic 
plate,  in  a  state  of  plane  strain. 
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<^1  = 


<^0 


_  2n 


l-2v 


(l-v)^  +  v| 
dr 


1  3up 


u, 


(1-v 


13ue 
r  30 


l-2v[ 

[1 3Ur  3Uq  Ufl ' 

r  d0  dr  r  J 


r 

r 


30 


+v- 


r  ; 

3Ur 

3r 


(8) 


on  91,  for  0  <  r  <  oo,  wherein  |j.  is  the  shear  modulus,  v  Poisson's  ratio"'^ ;  and  cohesive 
stress  boundary  conditions  on  0  =  0  together  with  any  of  the  admissible  set  of  boundary 
conditions  given  in  Table  1  on  0  =  <|)  for  0  <  r  <  «>.  Specifically,  we  are  interested  in  the 
local  behavior  of  the  fields  complying  with  the  foregoing  in  the  vicinity  of  the  wedge 
vertex,  O. 


Table  1.  Homogeneous  boundary  conditions  for  plate  edges 


No. 

Conditions 

Description 

(i) 

ae=kue,  T:re=k'Ui. 

Cohesive  stress-separation  laws 

(ii) 

Oq  =  0,  Xre=0 

Stress  free 

(iii) 

Ur  =  0,  ue=0 

Clamped/rigid  adhesive  contact 

(iv) 

ue=0,  Tre=0 

Symmetry /rigid  lubricated  contact 

(v) 

Ur  =  0,  Oq  =  0 

Antisymmetry 

For  plane  stress,  replace  d  by  v  /  (1  +  •u). 
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Several  comments  are  appropriate  at  this  point.  We  begin  this  commentary  by 
considering  the  general  nature  of  the  preceding  formulation,  then  focus  the  discussion 
on  each  of  tiie  boundary  conditions  in  Table  1  in  turn. 

First  with  respect  to  the  selection  of  plane  strain  (or  plane  stress)  field  equations, 
we  remark  that  results  foimd  can  be  expected  to  apply  to  axisymmetric  configvurations, 
and  even  three-dimensional  ones  provided  the  geometry  varies  smoothly  enough  in  tiie 
out-of-plane  direction  (Aksentian  [8]).  Second  re  the  choice  of  homogeneous  boimdary 
conditions,  we  observe  that  typically  singular  behavior  is  controlled  by  these  conditions 
rather  than  their  inhomogeneous  counterparts-^  Third  regarding  the  absence  of 
regularity  requirements  at  infinity  on  the  open  wedge  SR,  we  note  that  this  renders 
fields  complying  with  our  formulation  nonunique.  Since  the  principal  attribute  of  these 
fields  is  the  potential  characterization  of  all  possible  responses  at  the  wedge  vertex,  such 
a  lack  of  uniqueness  is  to  be  desired  rather  than  regulated  against. 

Turning  to  the  specific  boundary  conditions  in  Table  1,  the  first  set  is  that 
associated  with  cohesive  stresses.  Thus  k,k^  in  (i)  are  stiffnesses  associated  with  relative 
displacements  between  material  on  the  two  sides  of  the  ray  on  which  the  condition  is 
applied.  These  "springs"  need  not  be  linear;  in  fact,  typically  k  =  k(ue),  k'  =  k'(Ur). 
The  second  and  third  conditions  are  the  traditional  "free"  and  "clamped"  cases  treated 
by  Williams  [2];  the  latter  of  these  two  also  admits  to  interpretation  as  the  homogeneous 
complement  to  indentation  by  a  rigid  punch  with  complete  adhesion.  These  three  sets 
of  boundary  conditions  realize  three  problems  within  the  class  imder  consideration:  (i) 
"  (i)/  (i)  ■  (ii)/  si'd  (i)  -  (iii).  In  the  instance  wherein  the  cohesive  stress  conditions  are 
taken  on  both  faces,  it  is  useful  to  distinguish  between  symmetric  and  antisymmetric 
response.  The  next  two  sets  of  boimdary  conditions,  (iv)  and  (v),  enable  one  to  do  this. 
The  symmetry  conditions  can  also  be  interpreted  as  the  homogeneous  complement  to 

^  Usually  there  are  no  singularities  which  stem  from  the  inhomogeneous  part  of  physically  sensible 
boundary  conditions.  However,  jump  discontinuities  from  one  wedge  face  to  the  other  in  or  ue  /  r 

can  produce  Inr  stress  singularities. 
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indentation  by  a  perfectly  smooth,  or  lubricated,  rigid  punch.  In  all,  therefore,  the 
independent  problems  within  the  class  of  present  concern  are:  (i)  -  (ii),  (i)  -  (iii),  (i)  -  (iv), 
and  (i)  -  (v). 

2.  Equivalent  boundary  conditions  for  cohesive  stress-separation  laws 

In  this  section,  we  consider  the  role  of  cohesive  stress-separation  laws  as 
boundary  conditions  to  see  the  degree  to  which  it  can  be  reduced  to  that  of  other,  more 
traditional,  boundary  conditions.  In  this  way  we  hope  to  reduce  the  analysis  of 
cohesive  stress  boundary  conditions  by  taking  advantage  of  known  asymptotic 
solutions. 

A  sketch  of  a  normal  cohesive  stress  versus  separation  law  stemming  from 
attraction/repulsion  between  atoms/molecules  is  shown  in  Fig.  4.  Therein  s  is  the 
separation,  Se  its  equilibrium  value.  As  the  slope  varies,  we  can  expect  k  of  (i).  Table  1, 
to  be  a  function  of  separation.  However,  initially  as  loading  commences,  a  constant 
stiffness  is  a  fair  approximation  (as  indicated  by  the  dashed  line  in  Fig.  4).  Thus  we  take 
k  to  be  constant  henceforth.  For  similar  reasons,  we  can  also  approximate  k'  as  being 
constant. 

To  simplify  matters  further,  we  temporarily  take  k'  =  0  and  focus  on  the  normal 
cohesive  stress.  That  is,  we  consider 

a9  =  kue,  i:re=0  at  0  =  0,  (9) 

for  0  <  r  <  oo,  where  k  is  now  constant.  This  simplification  is  actually  an  appropriate 
one  if  cohesive  stresses  act  on  a  line  of  symmetry,  because  then  there  is  no  relative  Ur 
and  tj-e  =  0..  Notice,  too,  that  on  such  a  line  Uq  need  not  be  zero  if  cohesive  stresses  act 
(e.g..  Fig.  2),  in  contrast  to  the  symmetry  conditions  in  (iv).  Table  1. 
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Prior  to  actually  complying  with  (9),  we  need  to  construct  stress  and 
displacement  fields  satisfying  the  governing  equations,  (7),  (8).  Following  Knein  [1], 
Williams  [2],  we  use  the  Airy  stress  function,  X,  to  this  end.  That  is,  we  take 


throughout  91,  with 

V^X  =  0,  (11) 

on  91,  being  the  biharmonic  operator.  While  X  is  biharmonic,  it  can  be  obtained  as 
a  combination  of  harmonic  functions,  Vi  and  ¥2,  via 

X  =  Vi+rV2-  (12) 

The  key  in  choosing  ¥1  and  ¥2  is  to  ensure  that  the  resulting  X  has  four  linearly 
independent  functions  of  0  which  share  a  common  r-dependence:  then  when  boundary 
conditions  are  ultimately  applied  on  6  =  0,  <[),  we  will  have  available  four  arbitrary 
coefficients,  or  constants,  to  satisfy  the  four  conditions  involved.  With  this  in  mind,  we 
start  by  selecting  the  classical  separable  solution  for  ¥1, 


Vi  =  r^'^^[Acos(X  + 1)6  +  B  sin(X  + 1)6] . 


(13) 
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In  (13),  A  and  B  are  two  of  the  required  constants,  and  the  exponent  X.+1  is  chosen 
merely  for  subsequent  convenience  in  expressing  eigenvalue  equations  (after  Williams 
[2]).  Next,  to  obtain  two  more  functions  of  0  which  share  as  a  coefficient,  we 
exchange  X  + 1  for  X  - 1  in  (13)  and  take 

¥2  =  cos(^  - 1)0  +  D  sin(X  - 1)0],  (14) 

where  C  and  D  are  two  further  constants.  Combining  Yi  and  ^2  in  accordance  with 
(12)  then  yields  the  desired  %  with  a  total  of  four  constants  with  the  same 
multiplier.  The  attendant  stresses  can  be  derived  directly  from  the  Airy  stress  function 
relations  (10).  The  corresponding  displacements  can  be  similarly  determined  with  the 
aid  of  a  further  auxiliary  harmonic  function  (see  Williams  [2]  for  details),  or  simply  by 

integrating  the  stress-displacement  relations  (8).  The  resulting  fields  have,  for  the 
stresses, 

[(>.  -  3)(a  cos(X  - 1)0  +  b  sin(X  - 1)0)  +  c  cos(X + 1)0  +  d  sin(X  + 1)0], 

Oq  =  Xr^~^  [(X  +  l)(a  cos(X  - 1)0  +  b  sm(X  - 1)0) + c  cos(X  + 1)0  +  d  sin(^  + 1)0],  (15) 

Xr0  =  ^r^~^  [(A.  -  l)(a  sin(A  - 1)0  -  b  cos(A  - 1)0)  +  c  sin(A  + 1)0  -  d  cos(A  + 1)0], 

in  which  a,  b,  c,  d  are  constants,  and,  for  the  displacements, 

-r^ 

cos(A  - 1)0  +  b  sin(A  - 1)0)  +  c  cos(A  + 1)0  +  d  sin(A  + 1)0], 

,  (16) 

I 

ue  =  —  [(X  +  K)(a  sin( A  - 1)0  -  b  cos(A  - 1)0) + c  sin(A  + 1)0  -  d  cos(A  + 1)0], 
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wherein  k  =  3-4v.‘*’ 

Now  returning  to  the  satisfaction  of  the  cohesive  stress-separation  law,  the 
second  of  (9)  modifies  ae,Ue  of  (15),  (16)  so  that 

Oq  =  Xr^”^  [(A,  +  l)(a  cos(X,  - 1)0  +  b  sin(X.  - 1)6) + c  cos(A,  + 1)0  -  (A,  -  l)b  sin(A,  + 1)0], 

(17) 
r^ 

ue  =  —  [(X  +  K)(a  sin(X  - 1)0  -  b  cos(X  - 1)0) + c  sin(X + 1)6  +  (X  -  l)b  cos(X  + 1)0]. 

Substituting  (17)  into  the  remaining  condition  in  (9)  it  becomes  apparent  that  the  Cq  and 
U0  of  (17)  caimot  interact  because  of  their  different  powers  of  r.  At  first  thought,  then,  it 
would  appear  that  Oe  and  Uq  must  both  be  independently  zero  in  order  to  satisfy  the 
first  of  (9),  since  this  is  a  homogeneous  condition.  Unfortunately,  such  a  double 
requirement  can  seldom  be  complied  with  since,  in  effect,  it  realizes  a  total  of  five 
boundary  conditions  and  we  still  have  only  four  constants  to  meet  them  with. 
Fortunately,  this  overly  restrictive  situation  can  be  alleviated  as  follows.  We  do  indeed 
let  the  quantity  with  the  lower  power  of  r  be  independently  zero.  Thus  we  let  Gq  of  (17) 
by  itself  be  zero.  If  we  then  increase  X  by  one  so  that  X  -»  X  + 1,  the  associated  stress 
can  interact  with  Ue  of  (17).  However,  the  X  +  1  field  has  a  Ue  leftover;  this 
displacement  in  turn  must  interact  with  a  Gq  for  X->X  +  2.  And  so  on.  That  is, 
replacing  X  by  X  +  n,  the  first  of  (9)  requires: 

cfe  I  =0, 

n=0 

(18) 

del  =  kue  I  ,  n  =  l,2,... . 
n  n-1 


For  plane  stress,  k  =  (3  -  v)  /  (1  +  v). 
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Straightforward  manipulation  then  yields  the  following  series  expansions  for  <Je,Ue: 

<Te=ir^™-*  a„  £„(e), 

n=0 

(19) 

ue=ir‘*"  b„  £„(6), 

n=0 


where  the  vectors  an,  bn,  fn  (6)  are  given  by 

a„  =  (X  +  n)((>.  +  n  +  l)an ,  (X  +  n + l)bn ,  K^b^-i  -  (X  +  n + l)an  ,-{X+n-  l)bn ), 
=-^(-(?^  +  n+K)bn,(>.  +  n+K)an,(>.  +  n-l)bn,KX-i-(>-  +  n+l)an),  (20) 
fn(6)  =  (cos(X  +  n  - 1)0,  sin(X  +  n  - 1)0,  cos(X  +  n  + 1)0,  sin(^  +  n  + 1)0), 

with 


-k(K  +  l) 
2p(X  +  n)' 


(21) 


and  the  understanding  that  b-i=0  while  the  other  an,  bn  are  simply  the  constants  for  X 
being  +  n,  n  =  0, 1, 2,... . 

Substituting  the  series  expansions  of  (19)  into  whatever  are  the  outstanding 
boundary  conditions  that  hold  on  the  other  wedge  face  at  0  =  <|)  results  in  a  series  of 
2x2  systems  of  equations  in  the  constants  an,  bn  (n  =  0, 1,  2,...).  The  first  of  these  is 
homogeneous.  It  requires  the  value  of  X  be  an  eigenvalue  of  the  problem  associated  with 
^  and  relates  bo  to  ao-  Subsequent  systems  are  inhomogeneous  and  serve  to 
determine  an,  bn  (n  =  1, 2,...).+ 

For  n  =  0,  (9)  reduces  to  (see  (18)) 

^  This  is  so  provided  the  determinants  of  the  attendant  coefficient  matrices  are  not  zero  -  we  demonstrate 
how  to  handle  such  systems  in  Section  4. 


16 


ae  =  Xr0=O.  (22) 

It  follows,  therefore,  that  the  cohesive  stress-separation  conditions  (9)  have  an  identical 
influence  on  the  eigenvalues  of  one  of  our  wedge  problems  as  do  the  stress-free 
conditions.  Moreover,  this  can  similarly  be  shown  to  be  true  for  the  dual  of  (9),  namely 
k  =  0,  k'  0.  Together,  then,  these  results  establish  that  the  boundary  conditions  for 
cohesive  stress-separation  laws  ((i)  of  Table  1)  are  completely  equivalent  to  the  stress-free 
boundary  conditions  ((ii)  of  Table  1)  in  so  far  as  the  eigenvalues  of  any  of  our  admissible 
wedge  problems  are  concerned.  This  means,  in  particular,  that  any  singular  nature 
present  in  a  wedge  with  cohesive  stress  boundary  conditions  is  the  same  as  titat  present 
with  these  conditions  exchanged  for  stress-free  boimdary  conditions. 

3.  Eigenvalue  equations 

Here  we  exploit  the  equivalence  of  cohesive  stress  boundary  conditions  with 
stress-free  conditions  to  establish  a  complete  set  of  eigenvalue  equations  for  cohesive 
stress  boimdary  conditions  in  combination  with  any  of  the  other  boundary  conditions  in 
Table  1.  We  then  examine  these  equations  for  possible  singular  behavior. 

As  mentioned  in  Section  2,  when  cohesive  stress  boundary  conditions  hold  on 
9  =  0,  the  further  boundary  conditions  on  6  =  <|)  generate  a  homogeneous  2x2  set  of 
linear  algebraic  equations  in  ao,  bo.  For  a  nontrivial  solution,  the  coefficient  matrix  of 
this  set  must  be  singular.  Setting  the  determinant  to  zero  then  provides  an  equation  in 
X  such  that  this  is  so,  i.e.,  the  eigenvalue  equation. 

Given  the  equivalence  with  respect  to  eigenvalue  equation  generation  of  cohesive 
stress  boundary  conditions  ((i)  of  Table  1)  with  stress-free  conditions  ((ii)  of  Table  1),  the 
following  distinct  pairs  of  boundary  conditions  for  the  angular  wedge  may  be 
identified.  First  we  have  these  equivalent  conditions  paired  with  themselves: 

(i)  or  (ii)  -  (i)  or  (ii). 


(23) 
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Asymptotic  character  found  for  this  pairing  is  the  same  as  results  from  both  its 
symmetric  plus  its  antisymmetric  contributors: 

(i)  or  (ii)  -  (iv)  &  (v).  (24) 

Completing  possible  pairings  with  (i)/ (ii),  we  have  the  outstanding  set  of  boundary 
conditions  of  Table  I,  namely 


(i)  or  (ii)  -  (iii).  (25) 

The  three  eigenvalue  equations  resulting  from  (24),  (25)  represent  a  complete  set  of  all 
possible  combinations  of  the  boundary  conditions  in  Table  1  with  cohesive  stress 
boimdary  conditions.  These  equations  are  set  out  in  Table  2  (p.  18). 

The  eigenvalue  equation  for  (23)  is  given  in  Williams  [2];  separating  it  into 
symmetric  and  antisymmetric  parts  yields  the  first  two  equations  of  Table  2.  In  these 
equations,  <])  continues  as  the  entire  vertex  angle  so  that  (iv),  (v)  have  actually  been 
applied  on  6  =  <[)  /  2,  the  wedge  bisector.  The  eigenvalue  equation  for  (25)  for  plane 
stress  is  given  in  Williams  [2]  as  well;  adaptation  to  the  case  of  plane  strain  for  inclusion 
in  Table  2  is  immediate. 
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Table  2.  Eigenvalue  equations  with  cohesive  stress  boundary  conditions 


No. 

Equation 

Applicable  boimdary  conditions, 
associated  vertex  angle 

I 

sinXcj)  = -X,sin(t) 

(i)-(iv),  (t)/2;(i)-(i)or(ii),  <|) 

n 

sinXxj)  =  X,sin<|) 

(i)-(v),  <|)/2;(i)-(i)or(ii),  <|) 

m 

(k  + 1)^  =  4(X?  sin^  (j) + Ksin^  Xxj)) 

(i)  -  (iii),  q 

In  checking  for  roots  of  the  eigenvalue  equations,  we  are  primarily  concerned 
with  values  of  X  that  are  less  than  one  since  these  can  lead  to  singularities  (see  (15)  or 
(19)).  For  X,  =  1,  it  is  possible  to  have  a  logarithmic  singularity,  so  we  include  this  value 
too  (see  Dempsey  and  Sinclair  [9]  for  the  full  necessary  conditions  for  this  possibility). 
We  restrict  the  search  to  X.  >  0  so  that  any  such  stress  singularities  are  guaranteed  to  be 
integrable,  viz.,  to  have  finite  forces.  On  occasion,  eigenvalues  can  be  complex:  then  the 
preceding  restrictions  apply  to  the  real  part.  Thus,  in  sum,  we  confine  attention  to 
eigenvalues  in  the  range 

0<ReX.<l.  (26) 

Eigenvalues  complying  with  (26)  for  the  three  eigenvalue  equations  of  Table  2  are 
presented  in  Fig.  5  as  a  function  of  wedge  vertex  angle. 

The  curves  in  Fig.  5  display  values  of  the  singularity  exponent,  1-X,  rather  than 
the  eigenvalue  itself,  X..  The  larger  the  value  of  this  exponent,  the  more  singular  are 
possible  vertex  stresses.  This  is  because 


Singularity  exponent,! 


Fig.  5.  Singularity  exponents  for  cohesive  stress  boundary  conditions 


wherein  a  is  any  stress  component,  Of/Cfe/  or  x^e  (see  (15)).  Apparent  in  Fig.  5,  is  that 
the  symmetric  response,  I  of  Table  2,  is  uniformly  more  singular  than  the  antis5mimetric, 
n,  for  (j)  <  271.  Hence  I  is  the  curve  shown  in  Williams  [2]  for  the  dominant  singular 
behavior  imder  free-free  conditions.  Here  it  is  also  the  dominant  singular  curve  for 
cohesive  stress-separation  conditions  in  concert  with  either  themselves  or  stress-free 
conditions.  For  cohesive  stress-separation  conditions  together  with  clamped  conditions, 
m,  the  curve  for  d  =  1  /  4  closely  follows  that  given  in  Williams  [2]  for  clamped-free 
conditions  (the  u  used  in  [2]  for  plane  stress  corresponds  to  u=  3/13  here  for  plane 
strain).  The  additional  case  of  0)=  1/2  is  included  for  these  conditions  since  it  is  more 
singular.  For  HI,  for  both  values  of  u,  there  exist  other  curves  which  represent  less 
singular  response  than  that  shown  in  Fig.  5  for  a  given  wedge  angle  <|),  yet  nonetheless 
potentially  singular  stresses. 

As  well  as  giving  the  strength  of  possible  stress  singularities,  the  eigenvalue 
equations  determine  when  singular  behavior  is  no  longer  possible.  Identifying 
boundary  conditions  as  in  Table  1,  these  nonsingular  ranges  of  vertex  angles  are  as 
follows.  For  (i)  with  (i)  or  (ii)  under  symmetric  loading,  i.e.,  from  I, 

0<<1)<7C.  (28) 

For  (i)  with  (i)  or  (ii)  imder  antisymmetric  loading,  i.e.,  from  n, 

0<<l)<<t>*,  tan(l)*  =  (t)*(0<(t)*<2;i).  (29) 


From  Abramowitz  and  Stegim  [10],  p.  224,  (j)*  =  257.45°.  For  (i)  with  (iii),  i.e.,  from  in. 
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0<(t><sin  ^Vk  +  1/2, 

with  the  understanding  that  the  inverse  sine  takes  its  principal  value.  For  u  =  1/4, 1/2, 
K  =  3  -  4u  takes  on  the  values  2, 1,  and  (30)  gives  <|)  ^  jt  /  3,  jc  /  4,  respectively,  as  shown 
in  Fig.  5. 

Given  the  equivalence  of  boundary  conditions  for  cohesive  stress-separation  laws 
with  stress-free  conditions,  the  above  ranges  are  the  same  as  for  free-free  and  clamped- 
free  wedges.  How,  then,  does  the  use  of  cohesive  stress  boundary  conditions  lower  the 
number  of  instances  of  singular  behavior?  There  are  two  ways.  For  free-free  wedges, 
one  can  mtroduce  cohesive  stress  conditions  on  an  interior  ray  within  a  given  wedge,  a 
ray  upon  which  it  is  reasonable  to  consider  the  possibility  of  failure:  thus  the  effective 
wedge  angle  is  reduced  and  more  likely  to  fall  within  the  nonsingular  ranges  of  (28), 
(29).  For  clamped-free  wedges,  one  may  be  able  to  replace  the  clamped  conditions  with 
cohesive  stress  conditions:  hence  the  range  of  nonsingular  wedge  angles  goes  from  (30) 

to  (28)  or  (29)  and  nonsingular  response  is  more  likely.  We  demonstrate  these  two  ways 
in  the  next  section. 

4.  Eigenfunctions  for  some  special  cases 

Here  we  consider  the  effects  of  introducing  cohesive  stress-separation  law;s  for 
some  specific  local  features.  The  particular  configurations  examined  are:  a  crack  tip 
imder  symmetric  (mode  I)  loading  (Fig.  6a),  a  reentrant  corner  subjected  to  symmetric 
loading  (Fig.  6b),  an  epoxy-steel  butt  joint  under  tension  (Fig.  6c),  and  a  junction  of  three 
phases  in  a  titanium  aluminide  microstructure  (Fig.  6d).  In  what  follows,  we  treat  each 
of  these  instances  in  turn,  and  compare  local  stresses  with  and  without  cohesive  laws. 

For  the  perfectly  sharp  crack  under  mode  I  loading,  symmetry  enables  attention 
to  be  confined  to  the  upper  half-plane.  For  this  half-plane,  classical  boundary 
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conditions  take  stress-free  conditions  on  the  crack  flank,  symmetry  conditions  ahead  of 
the  crack.  That  is,  for  the  coordinate  system  of  Fig.  6a, 


=  '^re  =  0  at  9  =  jc, 

U0=O,  Tr9  =  0  at  0  =  0, 


(31) 


for  0  <  r  <  oo.  Using  the  equivalence  with  respect  to  eigenvalue  equations  of  stress-free 
boimdary  conditions  with  cohesive,  the  associated  eigenvalue  equation  may  be  taken 
from  Table  2  as  I  with  <|)  =  27c,  and  is  : 

sin  2nX  =  0.  ^^2) 

Consequently  the  dominant  singularity  admitted  (recall  (26))  has  X  =  1  /  2  and 

a  =  0(l/'\/r)  as  r-^0,  (33^ 

wherein  c  continues  as  any  stress  component.^"  This  is  the  classical,  inverse-square- 
root,  stress  singularity  of  present  day  fracture  mechanics. 

With  Barenblatt's  approach  [4],  cohesive  stresses  are  introduced  on  the  crack 
flank  while  symmetry  conditions  are  maintained  ahead  of  the  crack.  Given  symmetry, 
these  boundary  conditions  are 


Oe-~2kue,  -1^0=0  at  0  =  jc, 

(34) 

ue=0,  1:^0  =  0  at  0  =  0, 

Details  of  these  local  fields  are  well  known  and  may  be  found,  for  example,  in  Rice  [7]  at  p.  216. 
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for  0  <  r  <  oo.  The  minus  sign  in  (34)  is  simply  an  outcome  of  applying  the  cohesive 
conditions  on  0  =  tc,  and  the  2  reflects  the  total  separation  of  Ihe  crack  faces  under 
symmetric  loading.  In  the  light  of  the  equivalence  of  the  cohesive  conditions  in  (34)  to 
the  stress-free  conditions  in  (31),  we  have  the  same  eigenvalue  equation  for  this 
configuration,  to  wit  (32),  and  same  singular  character,  to  wit  (33).  Now,  though,  there 
are  two  instigators  of  singular  stress  fields:  one  being  far-held  loading,  the  other  being 
the  cohesive  stresses  on  the  crack  flank.  Since  these  share  the  same  singular  character, 
they  can  be  adjusted  to  cancel  one  another,  as  in  Bareriblatt  [4]. 

Here,  instead,  we  introduce  cohesive  laws  both  behind  and  ahead  of  the  crack  tip 
(as  in  Fig.  2).  That  is,  our  local  boimdary  conditions  for  symmetric  loading  are 

ae=-2kue,  Xre=0  at  0  =  jc, 

(35) 

Gq  =  2kue,  tje  =  0  at  0  =  0, 

for  0  <  r  <  oo.  The  eigenvalue  equation  for  these  conditions  is  still  I  of  Table  2,  but  now 
with  <t)  =  7C,  i.e., 

sin  Xn  =  0.  (36) 

This  is  an  example  of  the  reduction  in  effective  wedge  angle  mentioned  at  the  close  of 
Section  3.  Eigenvalues  with  integrable  stresses  as  a  result  become 

X  =  0,1,2,...  (37) 

None  of  the  associated  stress  fields  are  singular  (the  X  =  0  field  corresponds  to  a  rigid 
body  displacement).  However,  we  do  have  the  situation  mentioned  earlier,  namely  that 
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successive  coefficient  matrices  for  assembling  eigenfunction  series  expansions  are 
themselves  singular.  Fortunately,  by  combining  rigid  body  displacement  fields  with 
simple  polynomial  solutions  for  the  stresses,  we  can  construct  eigenftmctions  which  are 
not  series.  These  are  most  succinctly  e>q)ressed  in  rectangular  coordinates  (Fig.  6a),  and 
take  the  following  forms:  for  the  stresses. 


Oy  =  2k| 


Vo+Vi^  +  V2Pj  +0(x^) 


= -2kv2  — 
^  a 


(38) 


as  X  — >  0,  y  /  a  « 1,  while  for  the  displacements 


v  = 


X 

Vq  +Vi  — 

a 


1+ 


(1  +  K)ky 
4p 


V2 


2  3  —  K  2 

X  +- - + 


^(3(l+K)x2+(3-K)y2) 


1  +  K 


+  0(xy^)-f  O(x^), 


(39) 


(3  -  K)kx  Vi 
u  =  -vo^^ - L - L 


4p 


y  +  ^((3  -  k)x2  +  (1  +  K)y2) 


2V2X 


y+^(3(3-K)x2  +  (l  +  K)y2) 


+  0(y^)+0(x^y). 


as  X  ->  0,  y  /  a  « 1.  In  (38),  (39),  2a  remains  the  crack  length,  vq,  V|,  V2  are  arbitrary 
displacements  in  the  y-direction,  and  the  shear  stress  associated  with  these  coefficients 
is  everywhere  zero.  This  is  not  the  case  in  general,  however.  For  example,  the 
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companion  shear  stress  for  Cy  =kv3(x/a)^  has  terms  of  ord  and  ord  as 

X  0.  +  The  above  fields  may  be  verified  directly  by  substituting  them  into  the  stress 
equations  of  equilibrium,  the  stress-displacement  relations,  and  the  boundary 

conditions  (35)  on  converting  to  rectangular  coordinates;  when  this  is  done,  the  fields  in 
(38),  (39),  together  with  x^y  =  0,  are  indeed  foimd  to  satisfy  all  these  requirements. 

They  can  also  be  augmented  by  the  fields  associated  with  Ox  =  constant,  Oy  =  Xjg,  =  0. 
With  or  without  these  last,  clearly  the  fields  in  (38),  (39)  are  not  singular. 

For  the  closely  related  configuration  of  a  90 "  reentrant  comer  subjected  to 
symmetric  loading  (Fig.  6b),  classical  boimdary  conditions  for  the  upper  half  of  the  plate 
have  stress-free  conditions  on  the  inclined  face  together  with  symmetry  conditions 
ahead  of  the  corner.  That  is 

<^0='^re  =  O  6  =  37c/4, 

(40) 

U0  =0,  Xre  =  0  at  6  =  0, 


for  0  <  r  <  oo.  Thus  the  eigenvalue  equation  is  (from  Table  2  as  I  with  ^  =  3%/  2), 


.  37rX,  - 

sin  —  =  K. 
2 


Solving  (41)  numerically  for  X  within  the  range  0-1  yields  X  =  0.544.  Hence  singular 
stresses  of  order 

o  =  0(r'^  '^^)  as  r-^O  (42) 

Herein  the  ord  notation  implies  that,  if  f(x)  =  ord(ax**)  as  x  0,  a  and  a  being  independent  of  x,  then 

in  fact  f(x)  =  kax®  as  x  0  with  the  multiplicative  constant  k  not  being  zero.  This  is  in  contrast  to  tfie  O 
notation  which  admits  the  possibility  of  k  being  zero. 
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are  possible  in  the  corner.  Thiis  the  singular  character  here  is  only  a  little  weaker  than 
that  for  the  crack. 

Now  introducing  cohesive  laws  between  the  upper  and  lower  halves  of  the  plate 
with  a  reentrant  corner,  the  boundary  conditions  become: 

<^e  =-'Cre  =;^(r  +  Ur-U0)  at  e  =  37c/4, 

(43) 

ae=2ku0,  at  0  =  0, 

for  0  <  r  <  oo.  These  conditions  take  the  vertical  traction  on  the  corner  faces  to  be  k  times 
the  total  separation  of  the  faces  due  to  their  angle  and  to  any  vertical  displacement. 
They  also  take  the  horizontal  traction  to  be  zero  since  there  is  no  relative  horizontal 
displacement  between  the  upper  and  lower  halves  of  the  plate  because  of  symmetry. 
From  Table  2,  the  eigenvalue  equations  for  these  conditions  are  I  and  n  with  ^  =  3%/  A, 
viz.. 


sin 


37cX 

4 


(44) 


Again  the  wedge  angle  is  effectively  reduced.  Real  eigenvalues  with  integrable  stresses 
then  are 

^  =  0,1,  (45) 

the  first  of  these  corresponding  to  a  rigid  body  displacement  field.  Complex 
eigenvalues  are  also  possible.  These  occur  in  pairs  as  complex  conjugates.  The  complex 
pair  with  the  smallest  positive  real  part  may  be  determined  numerically  and  is 
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X  =  l  +  iT\,  ^  =  1.885,  Ti  =  ±0.361. 


(46) 


The  next  complex  root  has  %> 3,  and  further  roots  have  still  greater  real  parts.  Thus 
none  of  the  eigenvalues  of  (44)  lead  to  singular  stress  fields.  The  series  expansions  for 
the  associated  eigenfunctions  can  be  constructed  in  the  manner  outlined  in  Section  2. 
After  some  algebra,  this  yields:  for  tibe  stresses. 


Or  =  2kvo  sin^  6 


2 


Yl+Z^-ViVose±r3+3:^±^ 

Lv  a'  a ;  V  a'  a 


cos  36 


+2^1 + ^  j(sine  +  sin  30) 


+  ord(r^) 


V2Re 

V2lm 


j  [{X  -  1)(X  -  3)cos(^  - 1)6  - (X.2  -  l)cos(>.  + 1)6 

+(V2-X^  +  1)((A,  -  3)sin(A,  - 1)6  -  (X,  —  l)sin(X,  + 1)0)]  +  ord(r^ 


Cq  =  2kvo  cos^  6 


’  2 


3|1+I0_vi^) 
a'  ay 


cos0-^3  +  3^+ ^  j  cos  30 


+2^1+^j(3sin0-sin30) 


+ord(r2) 


(47) 


+  1^((X  + 1)  sin(X  - 1)0  -  (X,  - 1)  sin(X  +  1)0)  +  ord(r  ^  )|, 
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Xre  =  2kvo  sine  cose 


+{^2-X^  +  lj(cos(X  +  l)e  -  cos(A,  -  l)e)j + ord(r^) 

as  r  -»  0, 0  ^  e  <  Sjt  /  4,  while  for  the  displacements 
u,  =Vosine  +  ^^^^[K-l-2cos2e] 

+ord(r2) + ord(v2r^)  +  ord(v^r^), 

ue  =vocose  +  r  — ?X^sin2e  +  -^ 

L(l  +  K)a'  a  J 

+ord(r2)  +  ord(v2r^)  +  ord(v'2r^). 


as  r  0, 0  <  e  ^  371  /  4.  In  (47),  (48),  a  is  now  the  depth  of  the  reentrant  corner,  a'  is  a 
further  normahzing  length  defined  by  a'  =  4p/k(l  +  K),  Vo,Vi,V2  continue  as  arbitrary 
displacements  supplemented  by  v^,  and  X  satisfies  (44)  with  the  minus  sign  with  its 
actual  value  being  either  of  the  complex  roots  given  in  (46).  These  fields  may  be 
verified  directly  by  substitution  into  the  governing  field  equations,  namely  (7),  (8) 
together  with  stress  compatibility,  V2(a,  +  ae)  =  0,  for  the  higher  order  stress  terms. 
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and  the  boundary  conditions  (43).  When  this  check  is  performed,  the  fields  in  (47),  (48) 
are  indeed  found  to  comply  with  aU  these  requirements.  In  the  interests  of  brevity,  we 
terminate  explicit  expressions  for  series  terms  in  these  fields  at  ord  (r)  and  simply  give 
the  orders  of  the  next  terms:  other  than  the  somewhat  lengthy  algebra  involved,  there 
would  not  appear  to  be  any  impediment  to  determining  these  higher  order  terms 
explicitly.  As  before,  clearly  the  fields  in  (47),  (48)  are  not  singulat. 

Turning  to  the  qpoxy-steel  interface  where  it  meets  the  outside  free  surface  of  the 
composite  tensile  specimen  (Fig.  6c),  in  view  of  the  relative  rigidity  of  steel  compared  to 
epoxy,  classical  boundary  conditions  for  the  elastic  response  of  the  epoxy  are  clamped- 
free: 


Uj.  =  U0  =  0  at  6  =  Ji  /  2, 

(49) 

for  0  <  r  <  oo.  If  the  composite  specimen  is  thick  in  the  out-of-plane  direction  in  Fig.  6c, 
the  appropriate  field  equations  are  those  of  plane  strain,  i.e.,  (7),  (8).  Even  if  the 
specimen  is  the  normal  cylindrical  one  used  in  practice,  these  field  equations  are  still 
appropriate  for  examining  local  behavior  (see  Aksentian  [8]).  Hence,  given  the 
equivalence  of  stress-free  conditions  with  cohesive,  the  pertinent  eigenvalue  equation 
can  be  drawn  from  Table  2  as  HI  with  (j)  =  ti  /  2.  That  is, 

(k  + 1)^  =  4^?i^  +  K  sin^  ^ 

For  epoxy,  a  representative  value  of  Poisson's  ratio  is  3/8  corresponding  to  k  — 3/2. 
Then  solving  (50)  numerically  for  X  within  the  range  0-1  furnishes  X  =  2/3,  and  the 
dominant  singularity  possible  has 
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a  =  )  as  r  0. 

This  is  a  weaker  singularity  than  for  either  of  the  previous  two  configurations,  but  quite 
a  strong  one  nonetheless. 

Now  instead  we  introduce  a  cohesive  law  between  the  epoxy  and  relatively  rigid 
steel.  The  resulting  boundary  conditions  are  taken  to  bet 

ae  =  -kue,  Ur=0  at  e  =  7c/2, 

(52) 

0e=Xre=O  at  0  =  0, 

for  0  <  r  <  oo.  In  (52),  we  have  only  incorporated  a  cohesive  stress-separation  law 
normal  to  the  interface  and  not  admitted  the  possibility  of  relative  tangential 
displacements  between  the  epoxy  and  the  steel.  This  choice  is  arguably  the  simplest 
model;  however,  one  could  quite  reasonably  also  introduce  a  cohesive  law  for 
tangential  displacements.  Given  the  equivalence  of  the  very  first  of  (52)  with  the 
condition  Cq  =  0,  these  conditions  in  their  entirety  realize  the  same  boimdary  conditions 

as  for  a  complete  half-plane  under  antisymmetric  loading.  From  Table  2,  the  eigenvalue 
equation  is  thus  II  with  (])  =  tc,  viz., 

sin>,^  =  0. 

Associated  eigenvalues  with  integrable  stresses  are 


>-=0,1,2,...  . 


(54) 
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The  first  of  these  continues  to  correspond  to  a  rigid  body  displacement  field.  Since  the 
other  eigenvalues  are  separated  by  integers,  we  have  the  situation  alluded  to  in  Section 
2  wherein  coefficient  matrices  for  assembling  further  terms  in  eigenfunction  expansions 
are  themselves  singular.  To  overcome  this  difficulty,  we  need  to  develop  auxiliary 
fields  to  those  of  (15),  (16);  this  is  done  by  successive  differentiation  of  the  expressions  in 
(15),  (16)  with  respect  to  X.  Details  of  the  fields  so  formed  are  set  out  in  the  Appendix. 
These  fields  enable  eigenfunction  expansions  to  be  constructed.  The  algebra  involved  is 
somewhat  lengthy,  but  straightforward.  The  end  result  as  far  as  the  stresses  go  is: 

Cj  =  kuo[jt(K  - 1)(1  +  cos  20)  +  (3  -  k)(20  +  sin20)] 

-fo(3-K)|  ^(2  +  k'(1-4k"0  +  202))-^-^  COS0 

_  IHo  (lo  -  k'(3  -  4k"0  +  60^ ))  -  ^ + 3^]cos30 

_  1HO(3;i_4kV'-40(1-k'))-— (k"-0)  sin0 

a 

_  iliO(37:  +  4K'K"-40(3  +  K'))-— (k"-30)  sin30 
.  |A  ^ 

_lnr[f 45iio  (i_  K')  -Hllcos0  +  ( 4^(3  +  k')  -3-^lcos30 

Ll  11  aj  [  \i  ^  aj 

44k'  ((k"  -  0)  sin  0  4-  (k"  -  30)  sin  30) 

J 

-2K'^(lnr)2(cos0  4- 3cos30)| 

4-ord(uor^(lnr)^)  4-  ord(uir^(lnr)^)  4-  ord(uir^  Inr), 
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00  =  kuo[7c(K  - 1)(1  -  cos 20)  +  (3  -  KX2e  -  sin20)] 
i  ^  f  ^ 


+kr(3  -  K  H 


(cos8-cos30) 


ku 

44k'k" — 2-0(3cose  -  COS30) 
Vi 

ku, 


+  4k'k"  -  40(3  +  k'))  -  (k"  -  0) 


ku 


+lnr 


(3jt + 4k'k"  -  40(3  +  k'))  -  ^(k"  -  30) 

^(3  +  k')  -  3-^\cos0  -  COS30) 

\i  ay 

ku 

+4k'-jj^(3(k"  -  0)sin0  -  (k"  -  30)sin30)] 

+6K'-^^(lnr)^(cos0  -  cos30)| 
+ord(uor^(lnr)^ ) + ord(uir^(lnr)^) + ord(U|r^  Inr), 


sin0 
sin  30 


Tj-e  =  -kuo[7c(K  -  l)sin20  +  (3  -  k)(1  -  cos20)] 


-kr(3-K) 


^(371  +  4k'(k"  -  0))  -  k"-^ 
I  ^  a 


(cos0-cos30) 


kuo  ui 


0(cos0-3cos30) 


kun 


^(6  -  K'  +  2k'0(2k"  -  0))  -  ^ + 

I  \i  2a  a 


sin0 


ku 


S-(10  -  3k'  +  2k'0(2k"  -  30))  -^  +  3-^ 
\i  ^  2a  a 


sin  30 


+lnr 


ku 

4k'— ^(k"(cos0  -  cos 30)  -  0(cos0  -  3cos30)) 


^  ku, 


4^(1 +k')-^ 
y  Vi  a 


sin0  + 


^4^(3 +  k')- 3^^ 


V- 


sin30 


-2k'  (In  r)^  (sin  0  -  3  sin  30)| 


+ord(uor^(lnr)^)  +  ord(uir^(lnr)^)+ ord(uir^  Inr), 


(55) 
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as  r  0,  0  <  0  <  7t  /  2.  Companion  displacements  are: 

Ur  =  Uo|cos0  +  ■^[7c(k  -  l)(2cos20  +  K  - 1)+ 2(3  -  K)(sin20  +  (k  - 1)0  )]| 
+ord(uor^(lnr)^  j  +  ord^u^r^  Inr)  +  ord(uir^ ), 

(56) 

U0  =  -Uo|sin0  +  ^[7t(K  - 1)  sin  20  +  (3  -  k)(1  -  cos  20  +  (1 + K)lnr)]| 

+-^  +  ord(uor^(lnr)^)  +  ord(uir^  Inr) + ord(uir^), 
a  '  ' 

asr->0,  O<0<7i/2.  In  (55),  (56),  uo,Ui,Ui  are  arbitrary  displacements  in  the  x- 
direction,  a  is  now  any  normalizing  length,  and  k,K',K"  ^re  given  by: 

k  =  k/7c(l+K),  k'  =  (3-k)(1+k)/8,  k"  =  7c(5 - k) / 2(3 - k).  (57) 


These  fields  may  be  verified  directly  by  substituting  them  into  the  governing  field 
equations  and  the  boundary  conditions  (52).  Performing  this  check  shows  that  the 
fields  in  (55),  (56)  in  fact  meet  all  of  these  requirements.  Clearly,  too,  the  stresses  of  (55) 
are  also  free  from  any  singularities. 

As  our  final  example  we  consider  the  stress  intensification  that  can  occur  within 
the  microstructure  of  titanium  aluminide  alloys.  For  the  TiAl-base  lamellar  alloy 
designated  as  K5  in  Larsen  et  al.  [11],  yunctions  of  phases  with  different  orientations  are 
not  imusual  (see  Figure  1  (b)  in  [11]).  Such  a  featiue  is  sketched  in  Fig.  6d.  Herein  the 
different  alignments  of  the  anisotropic  phases  effectively  presents  different  moduli  to 
the  interfaces  involved  when  response  is  in  the  elastic  regime.  In  Fig.  6d,  we  have 
simply  designated  these  different  moduli  as  M.i,Ui  (i  =  1,  2,  3)  for  each  of  the  three 
phases  included  in  our  example.  In  reality  the  situation  is  more  complex  than  this,  but 


35 


this  simplified  approach  serves  to  demonstrate  the  differences  between  response  with 
and  without  cohesive  laws. 

Classical  boundary  conditions  for  three  distinct  elastic  materials  sharing  three 
interfaces  assume  continuity  of  stresses  and  displacements  at  an  interface,  i.e.,  that  the 
materials  are  perfectly  bonded  to  one  another.  For  the  geometry  and  coordinates  of  Fig. 
6d,  these  conditions  may  be  expressed  as: 

I  =<^e  I  /  I  =tre  I  , 

e=o+  e=0"  e=o+  e=0" 

.  . 

Ur  I  =Ur  I  ,  ue  I  =ue  |  , 

e=0"  0=0'^  0=0” 

for  0  <  r  <  oo,  with  compaiuon  matching  conditions  on  6  =  (l)i,(l)i  +  <1)2.  In  (58),  6  =  0^ 
denotes  the  limit  9  0, 0  >  0,  while  0  =  0"  denotes  the  limit  0  ^  0,0  <  0.  Unfortunately 

no  full  asymptotic  analysis  of  such  a  trimaterial  jimction  is  available  in  the  literature. 
However,  by  considering  some  limiting  cases,  we  can  approximately  gauge  the  range  of 
possible  singular  behavior.  To  this  end,  set  Pi  =  |i2/  Vj  =  1)2  =  1)3,  then  consider  the 
two  limits 

Hs  /  Pi  0/  Ps  /  Pi  ^  (59) 

For  the  first,  P3  vanishing  lets  the  edges  of  the  combined  pj  -p2  wedge  approach  being 
"stress-free".  Thus  from  Williams  [2],  if  (j)i  +  (j)2  >  tc,  we  have  possible  singular  response 
of  the  form 


a  =  0(l/r^  ^)  as  r^O,  1/2<X<1. 


(60) 
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For  the  second,  ji3  becoming  absolutely  rigid  renders  the  edges  of  the  combined  |i|  —  \i2 
wedge  "clamped",  and  again  from  Williams  [2],  if  <|)i  +  <1)2  >  upper  value  of  <])  in  (30),  we 
have  singular  behavior  of  the  type  indicated  in  (60).  Note,  though,  that  the  precise 
values  of  X  entailed,  while  falling  within  the  range  given  in  (60),  are  not  generally  equal 
for  these  two  cases.  In  between  the  1I3  of  (59),  we  have  the  instance  of  IJ.3  =  pj  for  which 
there  is  no  singularity.  Hence,  in  all,  we  can  expect  singular  response  wititi  classical 
boundary  conditions  for  the  trimaterial  junction  depicted  in  Fig.  6d  to  vary  from  being 
as  strong  as  the  inverse-square-root  singularity  of  the  classical  crack,  to  being 
nonsingular,  depending  upon  material  properties. 

Now  we  introduce  cohesive  laws  on  the  interface.  That  is,  we  take 

cre=kue,  Xr0=k'Ur  on  0  =  0,  (61) 

for  0<r  <«,  with  analogous  conditions  on  0i  =<j)i,  <I)i+<j)2-  Herein  U0,Ur  are  total 
displacements  measured  across  the  interface.  Thus,  in  the  light  of  our  earlier  result  that 
these  conditions  are  equivalent  to  stress-free  ones,  we  in  essence  have  three  stress-free 
wedges  so  far  as  singularities  go.  Provided  each  of  these  "wedges"  has  a  proud  comer, 
i.e.,  (l)j  <  7C,  i  =  1, 2, 3,  then  the  response  is  nonsingular  -  recall  (28),  (29).+  Moreover,  this 
result  holds  irrespective  of  the  values  of  material  constants. 

Concluding  remarks 

The  examples  of  Section  4  represent  four  demonstrations  of  configurations  which 
promote  singular  stresses  when  treated  with  classical  bormdary  conditions,  yet  which 
are  nonsingular  when  cohesive  stress-separation  laws  are  introduced.  Other  examples 
exist.  Accordingly,  since  singular  stresses  are  rubbish  physically,  the  introduction  of 

+  In  the  event  that  one  or  more  of  the  "wedges"  has  a  reentrant  comer,  the  configuration  can  still  be 
rendered  singularity  free  by  tf\e  introduction  of  cohesive  conditions  within  such  wedges  on  rays  upon 
which  it  is  physically  appropriate  to  entertain  failure  (cf.,  the  earlier  reentrant  comer  example). 


37 


cohesive  boundary  conditions  would  appear  to  offer  a  means  of  significantly  improving 
modelling  in  solid  mechanics. 

The  imderlying  reasons  for  this  alleviation  of  singular  behavior  may  be  described 
as  follows.  Cohesive  stress-separation  laws  are  fundamental  to  continuum  mechanics 
in  that  they  represent  the  solid  state  physics  at  the  atomic/molecular  level  which 
ultimately  leads  to  the  constitutive  relations  for  the  continuum.  As  a  consequence, 
when  boundaries  in  the  continuum  approach  one  another  to  within  the  distances  over 
which  the  stress-separation  laws  are  active,  consistent  modelling  requires  the 
introduction  of  cohesive  stress-separation  laws  in  the  boimdary  conditions.  Such  is  the 
case  for  the  crack  flanks  of  the  mathematically  sharp  crack.  Not  to  insert  cohesive 
conditioirs  on  these  crack  flanks  is,  in  effect,  setting  the  stiffness  in  the  cohesive  law  to 
zero  there.  This  realizes  an  abrupt  discontinuity  in  the  effective  stiffness  from  that 
taken  ahead  of  the  crack,  and  it  is  this  nonphysical  discontinuity  which  in  turn  produces 
singular  stresses.  However,  merely  introducing  cohesive  conditions  on  the  crack  flanks 
and  taking,  for  example,  symmetry  conditions  ahead  of  the  crack  (in  essence  as  in 
Barenblatt  [4])  also  represents  the  introduction  of  a  discontinuity  in  cohesive  law 
stiffness.  This  is  because,  when  the  transverse  displacement  ahead  of  the  crack  is  set  to 
zero  in  the  symmetry  conditions,  the  stiffness  is  effectively  being  taken  as  infinite,  in 
contrast  to  its  finite  value  on  the  crack  flanks.  Again  such  a  nonphysical  discontinuity 
produces  singular  response.  Only  when  a  cohesive  law  is  introduced  consistently  along 
the  crack  plane  (as  in  Fig.  2)  are  such  discontinuities  avoided  and  nonsingular  stresses 
result.  The  same  sort  of  argument  is  directly  applicable  to  the  reentrant  corner  example 
treated  herein. 

In  general,  the  removal  of  singular  stresses  stems  from  the  consistent  employment 
of  cohesive  stress-separation  laws  as  boundary  conditions.  This  consistency  must  exist 
between  the  cohesive  stress-separation  law  and  the  constitutive  relations  of  the 
continuum,  as  well  as  within  the  cohesive  boundary  conditions  themselves. 
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Appendix 

Herein  we  furnish  details  of  the  auxiliary  fields  needed  for  the  construction  of 
eigenfunctions  for  the  epoxy-steel  butt  joint.  The  first  of  these  stress  and  displacement 
fields  follow  from  differentiation  of  (15)  and  (16),  respectively,  once  with  respect  to  X. 
This  gives,  for  the  stresses, 

Off  =  -r {A,  In  r[( A,  -  3)(a  cos(A.  - 1)0  +  b  sin( A  - 1)0)  +  c  cos( A  + 1)0  +  d  sin(A  + 1)0] 
+A0[(A,  -  3)(-a  sin(A,  - 1)0  +  b  cos( A  - 1)0)  -  c  sin(A.  + 1)0  +  d  cos(A.  + 1)0] 

+(2A,  -  3)(a  cos( A  - 1)0  +  b  sin(A  - 1)0)  +  c  cos(A  + 1)0  +  d  sin(A  + 1)0], 

Oe  =  r^“^{Alnr[(A  +  l)(acos(A  - 1)0 + bsin(A  - 1)0)  +  c  cos(A  + 1)0  +  d  sin(A  + 1)0] 
+A0[(A  + 1)(— a  sin(A  - 1)0  +  b  cos(A  - 1)0)  -  c  sin(A  + 1)0  +  d  cos(A  + 1)0] 

+{2X  +  l)(a  cos( A — 1)0  +  b  sin( A  - 1)0)  +  c  cos(  A  + 1)0  +  d  sin(  A  +  1)0], 

Xre  =  r^"^  {A  lnr[(A  -  l)(a  sin(A  - 1)0  -  b  cos(A  - 1)0)  +  c  sin(A  + 1)0  -  d  cos(A  + 1)0] 
-i-A0[( A  -  l)(a  cos(A  - 1)0  +  b  sin(A  - 1)0)  +  c  cos(A  + 1)0  +  d  sin(A  + 1)0] 

+(2A  -  l)(a  sin( A  — 1)0  -  b  cos( A  - 1)0)  +  c  sin( A  + 1)0  -  d  cos(A + 1)0], 

in  which  a,b,c,d  are  constants,  the  bar  atop  them  serving  to  denote  that  they  need  not 
be  the  same  as  their  antecedents,  a,  b,  c,  d.  The  corresponding  displacements  are 

-r^  _  - 

Uf  =  -^{lnr[(A  -  K)(a  cos(A  - 1)0  +  b  sin(A  - 1)0)  +  c  cos(A  + 1)0  +  d  sin(A  + 1)0] 

+0[(A  -  K)(-a  sin( A  - 1)0  +  b  cos( A  - 1)0)  -  c  sin( A  + 1)0  +  d  cos( A  + 1)0] 

+a  cos(A  - 1)0  +  b  sin(A  - 1)0], 
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-  - 

U0  =  — {lnr[(?i  +  K)(a  sin(A,  - 1)0  -  b  cos(^  - 1)0)  +  c  sin(X,  + 1)0  -  d  cos(X  + 1)0] 

+0[(X  +  K)(a  cos(A,  - 1)0  +  b  sm(X.  - 1)0)  +  c  cos(^  + 1)0  +  d  sin(X,  + 1)0] 

+asin(A,  - 1)0  -  bcos(X,  - 1)0}. 

That  these  fields  in  fact  satisfy  the  governing  equations,  (7),  (8),  may  be  verified  by 
direct  substitution.  Satisfaction  of  these  field  equations,  however,  has  to  be  the  case 
since  (15),  (16)  comply  with  (7),  (8)  for  all  X,  and  (7),  (8)  do  not  themselves  involve  X.  It 
is  of  note  that  these  fields  contain  terms  that  are  distinct  from  the  original  ones  in 
Williams  [2],  and  which  are  not  included  in  the  classical  solution  of  Michell  (see  e.g. 
Timoshenko  and  Goodier  [12],  p.  133,  or  even  the  extended  version  of  Michell's  solution 
given  in  Little  [13],  on  pp.  166, 167). 

Further  differentiation  gives  the  following  stresses 

Or  =  -r^-HMlnr)2  [{X  -  3)(acos(^  - 1)0 + b  sin(>.  - 1)0)  +  ccos(X  + 1)0  +  dsin()i  + 1)0] 
+2A,0  In  r  [(X  -  3)(-a  sin(X,  - 1)0  +  b  cos(X  - 1)0)  -  c  sin(X,  + 1)0  +  d  cos(X  + 1)0] 
+2  In  r [(2 A,  -  3)(a  cos(A,  - 1)0 + b  sin( A  - 1)0)  +  c  cos(A,  + 1)0  +  d  sin(?i  + 1)0] 
+20[(2A  -  3)(-a  sin(A  - 1)0 + bcos(A,  - 1)0)  -  csin(A,  + 1)0  +  dcos(A,  + 1)0] 
-A0^[(A.  -  3)(acos(A  - 1)0  +  bsin(A,  - 1)0)  +  ccos(A,  + 1)0  +  dsin(A,  + 1)0] 

+2(a  cos(  A  - 1)0  +  b  sin( A  - 1)0)}, 

Oe  =  r  {A(ln  r)^  [(A + l)(a  cos(X  - 1)0 + b  sin(X  - 1)0)  +  c  cos(A  + 1)0 + d  sin(A  + 1)0] 

+2 A0  In  r [( A + l)(-a  sin(A,  - 1)0 + b  cos(  A  - 1)0)  -  c  sin(A, + 1)0  +  d  cos(A, + 1)0] 
+21nr[(2A + l)(acos(X  - 1)0 + bsin(X  - 1)0) + ccos(A,  + 1)0  +  dsin(A. + 1)0] 
+20[(2X  +  l)(-asin(A  - 1)0 + bcos(A,  - 1)0)  -  csin(A. + 1)0  +  dcos(A  + 1)0] 
-A.6^  [(A, + l)(a  cos(A  - 1)0 + b  sin(A  - 1)0) + c  cos(A  + 1)0  +  d  sin(A. + 1)0] 

+2(a  cos(  A  -1)0  +  6  sin(A,  - 1)0)}, 
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t,e  =  ^{X(lnr)2[(X - lXasm(X - 1)9  - bcos(X - 1)0)  +  csin(X  + 1)0  - dcos(X  + 1)0] 

+2X0  lnr[(A.  -  IXa  cos(X  - 1)0 + b  sin(X  - 1)0) + ccos(X + 1)0  +  dsin(X + 1)0] 
+21nr[(2X  -  l)(asin(X  - 1)0  -  bcos(X  - 1)0)  +  csin(X  + 1)0  -  dcos(X  + 1)0] 
+20[(2X  -  l)(a  cos(X  - 1)0  +  b  sin(X  - 1)0) + c  cos(X  + 1)0  +  d  sin(X  + 1)0)] 
-X0^  [(X  -  IXa  sin(X  - 1)0  -  b  cos(X  - 1)0)  +  c  sin(X  + 1)0  -  d  cos(X  + 1)0] 
+2(a  sin(X  - 1)9  -  b  cos(X  - 1)0)}, 


A  ^ 

in  which  a,b,c,d  are  also  independent  constants.  The  associated  displacements  are 
-r^ 

Ur  = — {(lnr)2  [(X  -  K)(a  cos(X  - 1)0  +  b  sin(X  - 1)0)  +  ccos(X + 1)0 + dsin(X + 1)9] 

+201nr[(X  -  K)(-asin(X  - 1)0 + b  cos(X  - 1)0)  -  c  sin(X + 1)0  +  dcos(X + 1)0] 
+2  Inr  [a  cos(X  - 1)0  +  b  sin(X  - 1)0]  +  20[-a  sin(X  - 1)0 + b  cos(X  - 1)0] 

-0^[(X  -  K)(acos(X  - 1)0 + b  sin(X  - 1)0) + c  cos(X  + 1)0 + dsin(X + 1)0]}, 

r^ 

ue  =  — {(lnr)2[(X  +  K)(asin(X  - 1)9  -  bcos(X  - 1)0)  +  csin(X  + 1)0  -  dcos(X  + 1)0] 

+201nr[(X  +  K)(acos(X  - 1)0  +  bsin(X  - 1)0)  +  ccos(X  + 1)0  +  dsin(X  + 1)0] 
+21nr[a  sin(X  - 1)0  -  b  cos(X  - 1)0] + 29[a  cos(X  - 1)0  +  bsin(X  - 1)0] 

-0^  [(X  +  K)(a  sin(X  - 1)9  -  b  cos(X  - 1)9)  +  c  sin(X  + 1)0  -  d  cos(X  + 1)0]}. 

Again  these  fields  may  be  verified  by  direct  substitution  in  (7),  (8).  The  foregoing 
auxiliary  fields  suffice  for  the  construction  of  the  eigenfunctions  up  to  the  terms 
explicitly  given  in  Section  4;  higher  order  terms  require  further  auxiliary  fields  which 
can  be  generated  by  yet  further  differentiation. 


